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New experimental data obtained mainly at the GSI/FRS 
facility allow to modify the empirical parametrization of frag- 
mentation cross sections, EPAX. It will be shown that minor 
modifications of the parameters lead to a much better re- 
production of measured cross sections. The most significant 
changes refer to the description of fragmentation yields close 
to the projectile and of the memory effect of neutron-deficient 
projectiles. 



I. INTRODUCTION 

The pioneering experiments of projectile fragmenta- 
tion at relativistic energies of 40 Ar and 48 Ca beams at 
the LBL Bevalac have demonstrated the potential 
of this method for the production of exotic nuclei. Based 
on these ideas, the SIS/FRS facility || at GSI has used 
also heavier projectiles like e.g. 58 Ni, 86 Kr, 129 Xe, and 
208 Pb to produce and study exotic nuclei For 
planning such experiments, when count-rate predictions 
are needed, analytical descriptions of fragmentation cross 
sections are useful. They are also useful in simulation 
programs for projectile-fragment separators (like e.g. IN- 
TENSITY §] or MOCADI §). Compared to physi- 
cal models of high-energy fragmentation reactions, which 
in general involve time-consuming Monte-Carlo calcula- 
tions, the virtue of an analytical formula lies in the short 
computing time and the possibility to calculate easily 
sub-microbarn cross sections that are beyond the reach 
of physical-model calculations. 

In 1990, Summerer et al. [jToj proposed a universal 
Empirical PArametrization of fragmentation CROSS sec- 
tions ("EPAX", Ref. which was based on and simi- 
lar to previous prescriptions by Rudstam fll|| and Silber- 
berg et al. |L2| . The parametrization was to a large extent 
based on multi-GeV proton-induced spallation cross sec- 
tions, since only scarce heavy-ion induced experimental 
data were available at that time. Meanwhile, more pre- 
cise data from relativistic heavy-ion-induced fragmenta- 
tion reactions together with recent results from projec- 
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tile fragmentation of heavy nuclei ( 197 Au and 208 Pb) on 
H2 targets |l^Jl^]) allow a more stringent comparison 
of proton- and heavy-ion induced isotope distributions. 
This comparison indicates that for heavy nuclei the two 
reactions lead to different isotopic distributions, which 
cannot be obtained from each other just by scaling fac- 
tors. This can be easily understood since heavy- ion in- 
duced reactions are expected to deposit more excitation 
energy in a nucleus than proton-induced reactions, mak- 
ing the final product distributions — after evaporation 
— broader and more neutron-deficient. Nevertheless, the 
data show that in both cases the isotopic yield distribu- 
tions can be well described by Gaussian-like analytical 
functions with parameters that vary smoothly as a func- 
tion of fragment mass [ ftpf . In the present paper, we will 
base the choice of these parameters exclusively on heavy- 
ion-induced reaction data. 

We will first review briefly the basic characteristics of 
the EPAX formula and then show which modifications 
are necessary to improve the accuracy with which the 
new experimental results can be reproduced. This will 
be followed by a brief comparison with similar attempts 
by other authors. 



II. THE EPAX FORMULA 

A. Basic characteristics 

The basic characteristics of the analytical description 
of high-energy fragmentation cross sections by the EPAX 
formula are the following |(J: 

• In the absence of systematic excitation-function 
measurements of heavy-ion induced fragmentation 
reactions, the formula is valid only for the so-called 
"limiting fragmentation" regime, i.e. for projec- 
tile energies where the fragmentation yields are no 
longer energy dependent, at least within the accu- 
racy of the formula (approximately within a factor 
of 2). This is certainly true for incident energies 
considerably above the Fermi energy in nuclei (« 
40 A MeV), in particular for the typical SIS ener- 
gies of 500 to 1000 A MeV. 

• The EPAX formula is meant to describe the frag- 
mentation of medium- to heavy-mass projectiles; 
nucleon-pickup cross sections are not included. No 
attempt is made to describe the fragmentation of 
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fissile nuclei. Therefore, the range of validity is lim- 
ited to projectiles from around argon to below the 
lead and bismuth isotopes. Predictions for produc- 
tion cross sections of fission products or of frag- 
ments below U where fission competition is signifi- 
cant require an elaborate description of the fission 
process, such as can be found e.g. in a recent pub- 
lication by Benlliure et al. |l5fl . 

• For fragments sufficiently far away from the pro- 
jectile (i.e. for mass losses larger than 15-20% of 
the projectile mass), the isotope distributions are 
largely independent of the original nucleus; their 
position, shape, and width depend only on the frag- 
ment mass number. This constitutes what has been 
termed the "residue corridor" and is related to the 
fact that the isotope distributions are mainly gov- 
erned by statistical evaporation from highly excited 
prefragments produced in collisions between rela- 
tivistic heavy ions. 

• For fragments that involve only a small mass loss 
from the projectile, the isotope distributions should 
be centered close to the projectile and their vari- 
ance should be small. Therefore, a smooth transi- 
tion is anticipated between the residue corridor and 
the projectile. The parametrization of this smooth 
transition constitutes the main task in designing 
the formula. 

In a first step, a parameter set has been searched for 
that describes the fragmentation yields from projectiles 
located close to the line of f3— stability. In a second step, a 
modification of the yield distributions due to the neutron 
or proton excess of projectiles located on the neutron- or 
proton- rich side of the line of f3— stability (the "memory 
effect" ) has been parametrized. 

B. The parameters of EPAX Version 1 

As explained in detail in Ref. Jl(| , the cross section (in 
barn) of a fragment with mass A and charge Z produced 
by projectile fragmentation from a projectile (A p , Z p ) im- 
pinging on a target (A t , Z t ) is written as 

a(A, Z)=Y A ■ a(Z prob - Z) 

= Y A ■ n -exp{-R- \Z prob -Z\ u <^) (1) 

The first term, Ya, represents the mass yield, i.e. the 
sum of the isobaric cross sections with fragment mass A. 
The second term describes the "charge dispersion", the 
distribution of elemental cross sections with a given mass 
around its maximum, Z prob . The shape of the charge 
dispersion is controlled by the width parameter, i?, and 
the exponent, U n (U p ), on the neutron- (proton-) rich side 
of the residue corridor. The factor n — y^R/ir simply 
serves to normalize the integral of the charge dispersion 
to unity. 



The mass-yield curve is taken to be an exponential as a 
function of A p — A. The slope of this exponential, P, is a 
function of the projectile mass. An overall scaling factor, 
S, accounts for the peripheral nature of fragmentation 
reactions and therefore depends on the circumference of 
the colliding nuclei: 

Y A = S ■ P ■ exp(~P ■ (A p - A)) (2) 

S = S 2 ■ {A^ 3 + Al /3 + S±) [barn] (3) 
In P = P 2 ■ A p + P 1 (4) 

The numerical values of the various constants can be 
found in Table | 

TABLE I. Constants used in EPAX Version 1 (Ref. @) 
and those used in Version 2 (this work). Note that in a few 
cases also the functional form to calculate a parameter has 
changed (see text). 
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The charge dispersion is characterized by the three 
parameters R, Z pro {,, and U. These three parameters 
are strongly correlated and difficult to obtain uniquely 
with a least-squares fitting technique. Note that the iso- 
bar distributions are not symmetric on the neutron- and 
proton-rich side, therefore U has two different values, U p 
and U n , on the proton- or neutron-rich side of the valley 
of /3-stability, respectively. In Ref. (hJ, the exponent U 
for the neutron-rich side of the isobar distribution was 
chosen as U n — 1.5, whereas the proton-rich side falls off 
like a Gaussian (U p — 2). The maxima of the charge dis- 
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persions, Z pro b, have been parametrized relative to the 
valley of /3-stability, 

Zprob = Zp + A (5) 

Zfj is approximated by the smooth function 

Z fj =,4/(1.98 + 0.0155 -A 2/3 ) (6) 

A is found to be a linear function of the fragment mass, 
A, for heavy fragments (A > A4), and is extrapolated 
quadratically to zero: 



A 



A 3 • A 2 if A < A 4 

A 2 • A + Ai if A > A 4 



(7) 



Similar to the parameter Z pro i> just discussed, the width 
parameter, R, is a function of fragment mass only, irre- 
spective of the projectile. In Ref. |[o| it was found that 
the experimental i?-values can be approximated by an 
exponential of the form 



In R = R 2 -A + R 1 



(8) 



The equations given above are sufficient to describe the 
"residue corridor", i.e. the yield distributions of projec- 
tiles located on the line of /3-stability if the fragment mass 
is far from the projectile mass. Close to the projectile, 
the following modifications have to be introduced [fTof: 



A = A-ll + di-iA/Ap-^) 2 } 
R = R-[l + n- (A/A p - r 2 ) 2 ] 



(9) 
(10) 



This serves to gradually reduce to zero the offset of Z pro b 
from the line of /3-stability and to decrease the width of 
the charge dispersion when A approaches the mass of the 
projectile, A p . 

A final correction applies if the projectile does not lie 
on the line of /3-stability. In this case, the A/Z ratio of 
the fragments will to some extent "remember" the A p /Z p 
ratio of the projectile ("memory effect"). For an analyti- 
cal description, it is simply assumed that the charge dis- 
persions are shifted by an amount A m which is a certain 
fraction of the distance of Z p from Zp p , the nuclear charge 
on the line of /3-stability for A p . Close to the projectile, 
this fraction is clearly close to unity (full memory effect), 
whereas it should gradually approach zero with increas- 
ing distance of the fragment mass A from the projectile 
(loss of memory) . The shape of the isobar distribution is 
assumed to be unchanged. Thus 

Zprob = Z, 3 + A + A m (11) 

where A m for neutron-rich projectiles is given by 

A m = [n x ■ (A/A t ) 2 + n 2 ■ (A/A t ) 4 ] ■ (Z p - Z Pp ) (12) 

The corresponding constants for proton-rich projectiles 
are termed p± and p 2 instead of n\ and n 2 . Numerical 
values for all constants are given in Table Q. 



With this parametrization, EPAX Version 1 was 
rather successful in describing the gross features of iso- 
tope distributions of high-energy projectile fragmenta- 
tion. This has been visualized e.g. in Refs. j|-(7|j where 
comparisons of EPAX with experimental data over many 
orders of magnitude in cross section can be found. As a 
particular example, we plot in Fig|l] experimental iso- 
tope distributions from 1 A GeV 208 Pb fragmentation (tJ 
in comparison with the EPAX Version 1 parametrization 
(dashed curves). 
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FIG. 1. Isotope distributions from 1 A GeV 208 Pb frag- 
mentation in a " ai Cu target in comparison with the old 
EPAX parametrization Version 1 (dashed curves) and the 
modified Version 2 (full curves). 

This set of data was also chosen to illustrate that the 
old EPAX version has problems to reproduce satisfac- 
torily the isotope distributions of very heavy fragments. 
As can be seen best for the low-Z isotope distributions, 
the dashed lines are centered too much on the neutron- 
rich side and exhibit too large maxima. Moreover, the 
cross sections for fragments close to the projectile (with 
masses A > 200) are predicted much too high. Minor 
discrepancies were also found for 129 Xe-fragment yields 
Q where the experimental distributions were found to 
be wider than the EPAX predictions. Serious deficien- 
cies were revealed when, compared to the EPAX pre- 
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diction, too large cross sections for neutron-deficient iso- 
topes from 58 Ni fragmentation and too small cross 
sections of neutron-deficient Sn isotopes from 112 Sn frag- 
mentation 1116(1 were measured. 



III. MODIFICATIONS OF THE EPAX FORMULA 

A. Fragmentation of projectiles close to /3-stability 

The availability of extensive cross-section measure- 
ments for projectile fragmentation of 40 Ar, 58 Ni, 86 Kr, 
129 Xe, and 208 Pb |17|j-0| down to /ib or nb cross sections 
allowed in a first step to find a better parametrization of 
the mass yields Y(A). In a second step, we adjusted the 
parameters describing the residue corridor, i.e. the width 
parameter, R, and the slope constants, U n and U p , of the 
quasi-Gaussian charge dispersion together with its cen- 
troid, Zproi,. In a third step, the correction factors for 
isotopic yields close to the projectile were modified. Fi- 
nally, modifications for projectiles outside the valley of 
/3-stability were redetermined. Numerical values for the 
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FIG. 2. Experimental mass yields, Y(A), for the reactions 

of 40 Ar+ 12 G Rrf 86 Kr+ 9 Be (b) Ref | ); 129 Xe+ 27 A1 

new constants are given in the fourth column of Table Q. (c, Ref. ||), and 208 Pb+ nc "Cu (d, Ref. Q) compared to the 

old EPAX parametrization Version 1 (dashed line) and to the 
new EPAX formula Version 2 (solid line). The open squares 
are the sum of the experimental cross sections for a given 
fragment mass A, whereas the open triangles represent the 
calculated cross sections using the new EPAX formula for 
isotopes missing in the experimental data for a given mass A. 
The full squares are the sum of the two values. 



1. Integrated mass and charge yields 



According to Eqs. (|2|-[lj), Y(A) is described by an expo- 
nential function depending on the fragment mass, A, the 
projectile mass, A p , and the target mass, A t . In EPAX 
Version 1, no additional correction close to the projec- 
tile has been performed. As shown by the dashed lines 
in Fig. this yields only a very rough agreement with 
experimental data. Note that the experimental mass 
yields had to be complemented by calculated cross sec- 
tions (from EPAX Version 2) where experimental data 
points were missing. For the 40 Ar, 86 Kr, and 129 Xe data, 
the additional calculated cross sections contribute only 
close to the projectile; in the 208 Pb case, however, the 
experimental data are less complete and more significant 
corrections were necessary. 

Fig. H shows that the slope somewhat further away 
from the projectile is reasonably well reproduced, how- 
ever, the absolute height as well as the slope close to the 
projectile mass needed to be modified. A much better 
overall normalisation was achieved by reducing the scal- 
ing factor S2 by a factor of 0.6, as demonstrated by the 
solid lines in Fig. |^. An additional correction introduced 
in the new EPAX version consists of an increase in the 
mass yield, Ya, (Eq. ||) close to the projectile according 
to 



Y A =Y A -[l + yi -{A/A p -y 2 ) 2 }. 



(13) 



for A/A p > j/2- These modifications give a satisfactory 
agreement with experimental data. In the case of 129 Xe, 
the slope far from the projectile mass deserves further 
attention. However, due to the lack of other complete 



experimental data, we were not able to significantly im- 
prove the mass-yield description. 

Another comparison of integrated cross sections from 
EPAX to experimental data can be performed for charge- 
changing cross sections. For this purpose, we use exper- 
imental data from Binns et al. |l9] . Fig. || shows a com- 
parison of the experimental data to the EPAX Versions 1 
and 2. As can be seen, nice agreement is achieved with 
the new EPAX formula, in particular close to the tar- 
get. For the 132 Xe data, the excess charge-changing cross 
sections for lower Z are most likely due to secondary re- 
actions in the targets which are unsufficiently accounted 
for. 



2. Parameters of the residue corridor 

In the upper part of Fig. |], we show, as a function 
of the fragment mass, A, the new position Z pro f>{A) of 
the "residue corridor", for 129 Xe and 208 Pb, fragmenta- 
tion. Data from 238 U fragmentation (Ref. pQ]) are also 
included in the Figure, they were not used to fit the pa- 
rameters Ai and can serve as a check. For clarity, we plot 
the offset, A, from Zp according to Eq. ||. The dashed 
lines are for EPAX Version 1; the full lines represent 
the new Version 2. Note that for better results close to 
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FIG. 3. Experimental charge-changing cross sections for 
132 Xe (a) and 197 Au ions (b) on 27 A1 targets in comparison 
with the old EPAX parametrization Version 1 (dashed line) 
and the new EPAX formula Version 2 (solid line) . The exper- 
imental data are from Binns et al. |lS[| . 

the projectile, we have added to Z pro \, according to Eq. |^ 
a tiny fragment-mass dependent correction amounting to 
0.002 -A. The lines sloping downward above A « 110 and 
A k 180 illustrate the return of Z pro h to /3-stability close 
to the 129 Xe and 208 Pb projectiles according to Eq. |[ 

In a similar way, the lower part of Fig. |^ visualizes the 
analytical description of the width parameter, R. While 
the functional form of the correction for A, Eq. |^, was 
kept identical with EPAX Version 1, the correction for 
R near the projectile is now written as 

R = R-[l + r 1 -A p -(A/A p -r 2 ) i ] (14) 

for A/A p > r%, which yields more narrow charge distri- 
butions for A close to A p . 

Only minor changes were introduced for the slope con- 
stants, U n and Up, of the neutron and proton- rich sides 
of the charge distributions, respectively. While U n was 
increased from 1.5 to 1.65 (constant for all A), U p is now 
slightly fragment-mass dependent with 

U p = Ui + U 2 ■ A + U 3 ■ A 2 (15) 

which approaches U p = 2.2 for large A, compared to a 
previous constant value of 2.0 (Gaussian curve). 



FIG. 4. a) Fragment-mass dependence of the offset. A, of 
Zprot from the line of /3-stability according to Eq. ra. The 
dashed curves denote the old EPAX parametrization, while 
the modified Version 2 is indicated by the full curves. The 
downward-curving lines above A m 110 and A « 180 illustrate 
the return of Z pro b to /3-stability close to the projectile (Eq. ^) 
for 129 Xe and 208 Pb, respectively, as a projectile. Data from 
fits of Eq. ^ to the results from Refs. ||[i],jL5[ are shown by 
the full squares, b) Fragment-mass dependence of the width 
parameter R according to Eq. |^. The upward-sloping dashed 
(full) curves indicate the shrinking of the widths near the 
129 Xe and 208 Pb projectile, respectively, for the old (new) 
version of EPAX according to Eq. [f^ (Eq. Q). Results from 
fits to the data from 129 Xe and 208 Pb fragmentation are shown 
by the full squares; those from 86 Kr fra gme ntation (triangles, 
Ref. Q]) have been approximated by Eq. (dash-dotted line). 

A close inspection of the width parameter, R, resulting 
from fits to the measured 86 Kr isotope distributions || 
shows that they are slightly but systematically decreased 
compared to Eq. || (see Fig. |J). This means that the 
isobar distributions are slightly wider for a neutron-rich 
projectile like 86 Kr than for a /3-stable one. Therefore, 
we tentatively introduce the following neutron-excess de- 
pendence of R: 

R' = R • [1 - 0.0833 ■ (Z(s p - Z p ) (16) 

where R is calculated according to Eq. ||. For neutron- 
deficient projectiles, R is not changed. 

For a typical example, the combined effect of the 
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FIG. 5. Isotope distributions from 800 A MeV 129 Xe frag- 
mentation in a 27 A1 target Q in comparison with the old 
EPAX parametrization Version 1 (dashed curve) and the 
modified Version 2 (full curve). 

changes described above is visualized in Fig. ||, which 
displays the isotope distributions from 800 A MeV 
129 Xe fragmentation measured over four orders of magni- 
tude In particular, the data measured near the max- 
imum are now much better reproduced. The stronger 
increase in R near the projectile mass, A p — 129, leads 
to a much faster decrease of the corresponding fragment 
yields at A « 124. Similar observations, though over a 
smaller vertical range, can be made for the fragment dis- 
tributions from 208 Pb fragmentation shown in Fig. [j]. It 
is interesting to note that even for 238 U fragmentation 
products in the mass range 100 < A < 130, produced 
in violent collisions and correspondingly long chains of 
statistical evaporation, rather good agreement between 
experimental data p0| and EPAX can be observed. 



B. Fragmentation of neutron- oder proton-rich 
projectiles 

In the previous EPAX version JXOj ] , the parametriza- 
tion of the "memory effect" was fitted, for the case of 
neutron-rich projectiles, to the results from 48 Ca frag- 
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FIG. 6. Examples visualizing the "memory effect" for neu- 
tron-rich projectiles impinging on 9 Be targets. Panels a) and 
b): S and F isotope distributions from 48 Ca fragmentation; c) 
and d): Se and Ni isotope distributions from 86 Kr fragmen- 
tation. The dashed curves are for EPAX Version 1, the full 
ones for the present Version 2. Experimental data are from 

Refs. as. 



mentation at 213 A MeV (see Eq. |l2]). For the case 
of 86 Kr, we were able to check the validity of this pre- 
scription in a different fragment-mass range J5| . We have 
kept this prescription for the current Version 2, since the 
agreement with experimental data has not changed sig- 
nificantly (see Fig. ^ for selected examples) . 

The situation is different, however, for the case of 
neutron-deficient projectiles. There, the previous ver- 
sion was unsatifactory from the beginning: Firstly, there 
were only very few radiochemical cross sections avail- 
able to derive a guess for the memory effect on the 
neutron-deficient side. Secondly, the resulting formula, 
A TO = [p2 ■ (A/A t ) 4 ] ■ (Z p — Zp p ), does not yield the lim- 
iting value A m /(Z p — Zp p ) = 1 for Z p = Zp p . These 
deficiencies became particularly obvious when Schneider 
et al. p6[ | measured only small cross sections for neutron- 
removal products from 112 Sn. This lead to a differ- 
ent functional form for the memory effect for neutron- 
deficient projectiles. The new parametrization reads 



exp [pi + p 2 ■ A/Ap] ■ (Zp - Z t 



/3p) 



(17) 



which is equivalent to a stronger but more rapidly 
decaying memory effect compared to the polynomial 
(Eq. |l2|) used previously. The consequences of the 
new parametrization of the memory effect for neutron- 
deficient projectiles are visualized in Fig. ^[ The small 
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dF _ d(log(a)) ^ -2 R 
~dZ ~ dZ ~ Zn(10) 



{Z — Z pro ij) 



(18) 



99 100 101 102 103 104 105 106 107 
fragment mass A 

FIG. 7. Experimental isotope distributions for Sn isotopes 
from neutron-deficient projectiles on 9 Be targets: a) from 
112 Sn fragmentation |lq| , b) from 124 Xe fragmentation 
plotted in comparison with the old EPAX parametrization 
(dashed curves) and the modified version (full curves). 



measured cross sections for 100 Sn [|T(| are now better 
reproduced, while the good agreement (over more than 
seven orders of magnitude!) with data measured in 124 Xe 
fragmentation pl[ | is maintained. 

We note in passing that Eq. [l?] is very similar to what 
has been suggested by von Egidy and Schmidt [^|| to 
reproduce production cross sections from antiproton an- 
nihilation. The overall agreement with the experimental 
fragmentation yields studied in our work is worse, how- 
ever, when we adopt their description of the memory ef- 
fect, in particular on the neutron-rich side. Moreover, we 
do not fit any set of parameters separately for a specific 
projectile as has been done in Ref. |^2| , but rather want 
to describe all systems with the same parameter set. 

The yield distributions from the neutron-deficient pro- 
jectile 58 Ni, studied in great detail by Blank et al. j|, 
deserve special attention. Not only do they represent a 
rare example of cross-section measurements down to the 
sub-nanobarn level, but also a case where severe discrep- 
ancies with EPAX Version 1 were observed: the yields 
of even-Z isotopes close to the proton drip line where 
found to be enhanced in experiment by factors of up to 
750 As a remedy for this deficiency, we have cho- 
sen to switch over from quasi-Gaussian to exponential 
slopes of the charge distribution above a certain gradient 
of the cross-section distribution. In order to achieve a 
steady transition between the two slopes, the transition 
point and the slope of the exponential have to be adjusted 
carefully. For this purpose, we calculate the derivative of 
the logarithm of the cross section (Eq. [j]): 



The transition point to the exponential slope, Z exp , 
can be calculated for the proton-rich side as a function 
of the fragment mass A according to 



dF 



Jexp{A) — Zp ro i,(A) t 



Zn(10) 
2 • R{A) 



(19) 



From Z exp on, the slope is exponential with the same 
gradient as Eq. [I] at this point. 

The gradient for which we switch to the exponential 
slope has to be parametrized as a function of the fragment 
mass A. As for the moment only the 58 Ni data exhibit 
this exponential trend (the measured cross sections in 
the 100 Sn region reach the same cross section level, but 
are further away from the proton drip line), we tried to 
adjust the function dF/dZ in such a way to reproduce 
the 58 Ni data without deteriorating the good agreement 
with the measured 124 Xe data. The function which fulfills 
reasonably well these criteria is the following: 



dF 

— = 1.2 + 0.647- (A/2) - 3 

(LZj 



(20) 



The result for the 58 Ni fragment-yield distributions is 
visualized in Fig. |^. The full line shows the new EPAX 
cross sections including the above modification of the 
slopes, whereas the omission of the exponential slope 
leads to the dotted cross-section curve. As a consequence 
of this modification, EPAX Version 2 predicts a produc- 
tion cross section e.g. for the doubly-magic nucleus 48 Ni 
of 4 • 10~ 13 b, which is within reach of an experiment. 

As stated earlier, the present adjustment for very 
proton-rich projectile fragments is only based on results 
from 58 Ni fragmentation Q. Therefore, caution is ad- 
visable when applying this parametrization especially to 
very light proton-rich fragmentation products. It would 
be interesting to compare the present parametrization to 
fragmention yields at the drip line from other proton- 
rich beams, e.g. from 36 Ar or 78 Kr. In the 100 Sn region, 
the present transition to an exponential slope has only 
a slight influence on the 100 Sn cross-section prediction 
from 112 Sn fragmentation (6.6 • 10 -12 b, to be compared 
to an experimental value of 1.2 • 10~ 12 b |l(|). It might 
be interesting to measure fragmentation cross sections 
for even more proton-rich nuclei in this mass region to 
compare to the modified EPAX formula. 



C. Overall quality of EPAX Version 2 

Besides the comparison between experimental data 
and EPAX predictions for individual elemental distri- 
butions, we tried also to visualize in a more global man- 
ner the overall quality of the new EPAX parametrization 
compared to the data and the improvements with respect 
to EPAX Version 1. 
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FIG. 8. Experimental isotope distributions for Ni to V iso- 
topes from 58 Ni+ 9 Be fragmentation Q in comparison with 
the old EPAX parametrization (dashed curves) and the new 
version (full curves) . The new EPAX formula without a tran- 
sition to the exponential slope is shown by the dotted curve. 
This effect is clearly visible below the 10 nb level. 




lo gio(<W n EPAX> lOflloKxptaEPAx) IOfllo( a ip.' EP«x) 



FIG. 9. Logarithm of the ratio between experimental frag- 
mentation cross sections and predictions of EPAX Version 1 
(dashed line) and the new EPAX Version 2 (solid line) for 
projectiles ranging from 40 Ar to 208 Pb. References to the 
experimental data are given in the text. 



For this purpose, we plot the cross section ratios 
&exp/&EPAX for all projectile-target combinations used 
in the present paper in Fig. ^. For most of the projectiles 
the logarithm of the cross-section ratio is centered around 
zero indicating an agreement on average between experi- 
ment and prediction. For the new EPAX Version 2, the 
sigma widths of these distributions vary between 0.5 (in 
the case of 40 Ar and 86 Kr) and 0.2 (in the case of 129 Xe). 
If we analyse all experimental data together, we obtain 
a sigma of 0.4. This demonstrates that the new EPAX 
formula can predict cross-sections in most cases within 
a factor of two. In almost all cases we observe a signifi- 
cantly smaller sigma for the new EPAX parametrization 
than for the old version. This is particularly striking in 
the cases of 58 Ni, 129 Xe, and 208 Pb. Only in the case 
of 48 Ca the agreement between experimental data and 
predictions deteriorates slightly, exhibiting now a shift 
to smaller predicted cross sections. It is likely that the 
memory effect for neutron-rich projectiles is the origin of 
this discrepancy. 

D. Comparison with other empirical 
parametrizations 

As has been mentioned in the Introduction and dis- 
cussed in detail in Ref. ]To| | , other empirical parametriza- 
tions have been proposed, mainly for proton-induced 
spallation reactions (e.g. Refs. [ p3p^ ]). Contrary to 
our formula, these parametrizations fit also the energy 
dependence of the cross sections. The former approach 
has been extended to describe heavy-ion induced spalla- 
tion reactions by scaling proton-induced cross sections by 
an energy-dependent factor [p5[ and has achieved good 
agreement with measured cross sections (in the 100 mb 



to 1 mb range) for fragments from medium-mass nuclei 
(i.e. for projectiles up to 56 Fe, Ref. P5|j ), 

For fragmentation cross sections from heavy-mass nu- 
clei, however, the Tsao et al. formula |^{| is less suc- 
cessful. This can be demonstrated e.g. by comparing in 
Fig. IhJ their prediction for the Pt isotope distribution in 




175 180 185 190 195 200 205 
fragment mass A 



FIG. 10. Platinum (Z=78) isotope distribution from 1 A 
GeV 208 Pb fragmentation in a nat Cu target jj in comparison 
with the new EPAX parametrization (full curve) and the Tsao 
et al. formula (Ref. |25j, dashed curve). 

the reaction 208 Pb + nat Cu to experimental data @. As 
we have mentioned in the Introduction, our physical un- 
derstanding of high-energy heavy-ion reactions suggests 
that proton-induced reactions produce more neutron-rich 
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isotope distributions of heavy elements than heavy-ion- 
induced reactions, therefore it is unlikely that even a bet- 
ter fit of the Tsao et al. formula to the data in Fig. [h] 
would apply also to the same isotope distribution formed 
in the p+ 208 Pb reaction. We believe that separate pa- 
rameter sets have to be fitted to the respective experi- 
mental data. For k1 A GeV protons impinging on heavy 
targets, the recent work of Refs. Jll|[l4|] provides for the 
first time a comprehensive data set that allows to extend 
previous work p3|]24[ | to heavier fragments. A complete 
parametrization of the bombarding-energy dependence, 
however, has to await more measurements at different 
energies for the heavy systems. 



E. EPAX predicitions for projectiles very far from 
stability 

The parameters of the new EPAX formula have been 
adjusted by fitting data from stable projectiles. As ex- 
plained above, one of the most difficult tasks of the 
present work was to reasonably parametrize the cross 
sections for fragments close to the projectile. If the pro- 
jectile is very far from the line of (3 stability, this task is 
even more difficult. 

The comparisons between EPAX predictions and ex- 
perimental data shown in the present paper indicate 
that the chosen parametrization is reasonable. A pos- 
sible check of its predictive power could be to compare 
the EPAX formula to other data from neutron-rich or 
neutron-deficient projectiles. However, only very few 
other high-quality data are available. 

As one check, we compare in Fig. [ll] the fragmenta- 
tion of a stable 24 Mg beam jl8| to a radioactive 28 Mg 
beam [^6| and of a stable 40 Ar beam Q] to a radioac- 
tive 43 Ar beam pq |. The elemental distributions for Mg, 
Na, and Ne isotopes shown in the left-hand column and 
those for S, P, and Si isotopes shown in the right-hand 
column nicely show the experimental memory effect and 
the high quality of the EPAX predictions. In a similar 
way we compared the results from fragmentation of 96 Ru 
and 96 Zr |^7| to EPAX predictions and found reasonable 
agreement. However, for projectile beams still further 
away from stability, the memory effect becomes more and 
more important. Therefore, we think that some caution 
should be applied when using EPAX predictions for pro- 
jectiles very far from stability like e.g. 132 Sn. Here, no 
experimental data at all exist to verify the EPAX pre- 
dictions. 



IV. SUMMARY AND CONCLUSIONS 

We have demonstrated that the quality with which 
EPAX reproduces measured high-energy fragmentation 
cross sections could be improved considerably by intro- 
ducing rather small modifications to the formula. This 
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FIG. 11. Experimental cross sections for the fragmentation 
of stable (open symbols) and radioactive beams (full symbols) 
compared to the new EPAX parametrization (dashed and 
full lines, respectively). On the left-hand side, results from 
24 Mg |l8| and 28 Mg |2q] fragmentation are shown, whereas 
on the right-hand side Ar 111 fragmentation is compared to 
43 Ar fed] fragmentation. 



makes EPAX a more reliable tool to predict production 
rates for secondary-beam experiments with medium- or 
heavy-mass exotic nuclei. 

Even after introducing these improvements, there are 
still some discrepancies with measured data which de- 
serve attention in future modifications of the EPAX for- 
mula. One aspect concerns the odd-even effects which 
were taken into account by von Egidy and Schmidt [p2[ . 
At present, we do not think that such a rather small 
modification (of the order of 20-30%) is necessary in view 
of the overall discrepancies of factors around 2 still ob- 
served, but it may become necessary once better pre- 
cision can be achieved. Another open question is the 
change in slope observed for very neutron-deficient 58 Ni 
fragmentation products. Here it would be desirable to 
measure systematically the formation cross sections for 
fragments close to the proton drip line for other medium- 
mass neutron-deficient projectiles, e.g. 36 Ar and 78 Kr. 

Up to now, we did not use medium-energy data (from 



9 



reactions at less then «100 A MeV) which are available 
from GANIL, RIKEN, or MSU. First of all, these data 
probably do not fulfill the condition of " limiting fragmen- 
tation" mentioned above. In fact, it is well known that at 
these energies nucleon-exchange reactions as well as pick- 
up reactions become increasingly important and thus al- 
ter in particular the cross sections close to the projectile. 
A striking example for this is the difference in experimen- 
tal cross sections for 100 Sn production from a 112 Sn beam 
of roughly a factor of 100 in recent GANIL and GSI ex- 
periments ]2l| , |28[ | . In addition, these intermediate-energy 
data are measured with separators having very low trans- 
missions (of the order of few percent) which are difficult 
to measure experimentally. It would be interesting, how- 
ever, to compare EPAX also to cross sections obtained 
at lower incident energies once high-quality data become 
available. 

A FORTRAN program of the EPAX formula as well 
as additional graphs comparing experimental data to 
EPAX pr edictions can be downloaded from the f ollowing 
address: ftp://ftpcenbg.in2p3.fr/pub/nex/epax . 
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